Quantum critical behavior of electrons at the edge of charge order 
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We consider quantum critical points (QCP) in which quantum fluctuations associated with charge 
rather than magnetic order induce unconventional metallic properties. Based on finite-T calculations 
on a two-dimensional extended Hubbard model we show how the coherence scale T* characteristic 
of Fermi liquid behavior of the homogeneous metal vanishes at the onset of charge order. A strong 
effective mass enhancement reminiscent of heavy fermion behavior indicates the possible destruction 
of quasiparticles at the QCP. Experimental probes on quarter-filled layered organic materials are 
proposed for unveiling the behavior of electrons across the quantum critical region. 
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Introduction. Quantum critical points occur at zero 
temperature second order phase transitions in which the 
strength of quantum fluctuations is controlled by an ex- 
ternal field such as pressure, magnetic field or chemical 
composition [l|. In recent years intensive studies have 
focused on itinerant electrons at the edge of magnetic 
order, being the heavy fermion materials J2]-|4J proto- 
typical examples. Much less explored but equally inter- 
esting are QCP arising from tuning the electrons close 
to a charge ordering instability. This situation is real- 
ized in the quarter-filled families of layered organic su- 
perconductors [i| of the a, /3" and 6'-(BEDT-TTF)2X 
types. Large electron effective mass enhancements and 
non-Fermi liquid metallicity at finite-T are observed in 
(MeDH-TTP)2AsF6 and K-(DHDA-TTP)2SbF6 above a 
critical pressure at which the charge order found at am- 
bient pressure melts [y, |7|. Such heavy fermion behavior 
may appear puzzling, considering the different 7r-orbitals 
of the organics as compared to the f-orbitals in the rare 
earths, but can find a natural explanation based on the 
universal properties of matter expected near a QCP. 

Charge ordering (CO) phenomena in quarter-filled lay- 
ered organic materials are observed in a wide variety 
of crystal structures, not limited to specific Fermi sur- 
face shapes or nesting. This indicates the importance 
of onsite and intersite Coulomb repulsion (8|, |9| between 
TT electrons as the driving force of CO, and in turn im- 
plies that electronic correlation effects similar to those 
found in half-filled systems are inevitably present. These 
should be considered together with the quantum critical 
fluctuations of the order parameter to understand the 
metallic properties in the neighborhood of the present 
Coulomb-driven transition. The latter can, in principle, 
differ from more standard charge density wave instabili- 
ties of the Fermi surface. 

In this Letter we analyze theoretically the possible ex- 
istence of a QCP at a CO transition driven by the quan- 
tum fluctuations associated with strong off-site Coulomb 



repulsion, in the absence of Fermi surface nesting. The 
influence of the T = singular quantum critical point 
on the finite temperature metallic properties is stud- 
ied based on finite-T Lanczos diagonalization [10|, |ll| 
of an extended Hubbard model on an anisotropic tri- 
angular lattice. At quarter filling {n = 1/2 hole per 
molecule), lattice frustration naturally leads to charge 
ordered metallic states with a single CO pattern. Our 
main result is the existence of a temperature scale T* 
in the dynamical and thermodynamic properties of the 
system, that is suppressed as the CO transition is ap- 
proached from the homogeneous metal. The electronic 
specific heat coefficient at low temperatures is found to 
be strongly enhanced close to the QCP in analogy with 
studies of quantum criticality in heavy fermions. Follow- 
ing this analogy, non-Fermi liquid behavior should oc- 
cur at finite-T in quarter-filled organic conductors at the 
edge of a CO instability and could be directly probed 
by several experimentally measurable quantities. Due to 
the universal character of the proposed quantum critical 
scenario, our results could also be relevant to broader 
classes of systems where Coulomb-driven CO occurs, not 
restricted to the particular ordering pattern studied here. 
Model. We focus on the extended Hubbard model: 
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on the anisotropic triangular lattice shown in Fig. [IJa), 
where tij = {tp,tc) are the transfer integrals between 
nearest neighboring molecules respectively along the di- 
agonal {p) and vertical (c) directions, Vij = {Vp,Vc) are 
the corresponding inter-molecular Coulomb interaction 
energies and U is the intra-molecular Coulomb repulsion. 
The model Eq. ([T|) has been studied via a variety of tech- 
niques due to its relevance to 0-type two-dimensional or- 
ganic conductors (see Ref . I4I for a recent review) . Here 



we follow the standard practice and neglect longer-range 
Coulomb interactions 12| as well as electron-lattice ef- 
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FIG. 1: (color online) (a) Lattice structure and microscopic 
parameters of the extended Hubbard model Eq. ((TJ. (b) 
Phase diagram obtained at T = from numerical diagonal- 
ization of Ns — 12 (full line) and 18 sites clusters (dashed). 



fects [l3| that are however essential to recover the var- 
ious ordering patterns realized in these materials. For 
the sake of simplicity we consider an isotropic inter-site 



repulsion Vc = Vp = V and set tc 



0, U 



t > 0. 



This choice is representative of the 6'-ET2X salts with 
X=CsCo(SCN)4, X=CsZn(SCN)4 and X=l3 where the 
molecular orbital overlap is strongly suppressed along 
the c direction [sl |l2[. These materials lie close to (on 
both sides of) the bandwidth controlled CO transition in 
Mori's phase diagram [8| and are therefore optimal candi- 
dates for the observation of an interplay between critical 
charge fluctuations and electronic correlation effects. 

Phase diagram. The phase diagram obtained at T = 
from the numerical diagonalization of the model Eq. 
([1]) on iVs = 12 and Ns — 18 site clusters is presented 
in Fig. [Ijb). The different phases can be identified by 
analyzing the behavior of the charge correlation func- 
tion NsC{q) = iV^i I]ii("»"i)e*'''^''- In the thermody- 
namic limit, this quantity diverges at a single wavevector 
Q 7^ at the onset of charge order. An accurate nu- 
merical determination of the phase boundaries relying 
on a proper finite-size scaling of the results is prohibitive 
for the fermionic system under study, due to the rapidly 
increasing size of the Hilbert space. We therefore iden- 
tify the r = ordering transition, Vco, as the locus of 
steepest variation of charge correlations upon varying the 
interaction parameters. An analogous procedure is used 
to determine the melting temperature of CO, Tco- In 
the physically relevant regime explored here, U/t < 20 
[14J the phase boundaries agree on the two cluster sizes. 

In the absence of nearest-neighbor repulsion, V — 0, 
the system remains in a homogeneous metallic phase up 
to arbitrary values of the local interaction U, as holes can 
effectively avoid each other at concentrations away from 
integer fillings. An instability towards a charge ordered 
state with 3-fold periodicity is realized instead upon in- 
creasing the inter-site interaction, V, as was previously 
obtained by different approaches l^, Il5l4l9| . The result- 
ing 3- fold ordering pattern is shown in Fig. [IJb). 



At low and moderate values of U/t < 5, down to U — 
0, the CO transition essentially follows the predictions 
of mean-field approaches I^, lla, [la]- A calculation in 
the random phase approximation (RPA) yields Vco =^ 
1.06t + U/6, which is shown as a dotted line in Fig. [Ijb). 
This law is correctly recovered by the numerical data at 
low U, but sizable deviations appear at as soon as U/t > 
10 due to the increasing effects of many-body electronic 
correlations. The boundary obtained numerically in this 
region is independent of the cluster size, and our value 
Vco/t = 3.5 at U/t = 10 is in good agreement with 
existing numerical results in larger systems 
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Before moving to the analysis of the correlated metal- 
lic phase at the edge of charge order, let us note that 
the charge correlation function also provides indications 
of a crossover taking place within the CO phase, sep- 
arating a conventional 3-fold state from a more exotic 
"pinball liquid" phase (l7|. The latter arises because 
at large U, mean-field like configurations where charge- 
poor molecules are completely depleted become energet- 
ically unfavorable, as these imply that each charge rich 
molecule should accomodate up to 3/2 holes on aver- 
age. To prevent double occupancy, part of the hole 
density necessarily spills out and decouples from the 
charge rich sublattice, resulting in a separate fiuid mov- 
ing freely in the remaining sublattice [17| . This partial 
ordering, occurring for V < U/3, corresponds to a value 
C(Q) — "-^/S, to be contrasted with the value C(Q) — ri^ 
obtained in the 3-fold state at large V . 

The correlated metal close to charge ordering. We 
start by analyzing the kinetic energy of the interact- 
ing system, a quantity that provides direct information 
on how the motion of the charge carriers is hindered 
by interactions, and can be evaluated with good accu- 
racy through finite-T Lanczos diagonalization. Its impor- 
tance in correlated systems has been recently recognized 
20, l2l| , and resides in the fact that it can in principle be 
accessed from optical absorption experiments, providing 
a quantitative measure of many-body correlation effects. 

The kinetic energy, K^ normalized to the non- 
interacting band value, Kq, is shown in Figs. [2][a) and 
(b), respectively for U/t = 5 and U/t = 15, for several 
values of the intersite repulsion across the CO transition. 
At U/t = 5 the kinetic energy at T = stays essentially 
unrenormalized, K/Kq > 0.9, upon increasing V all the 
way up to the CO transition occurring at Vco = 2.33i, 
as expected in a weakly correlated Fermi liquid. It then 
suddenly drops to a value K/Kq ~ 0.6 upon entering 
the charge ordered phase. This residual value is ascribed 
to local (incoherent) hopping processes in the charge or- 
dered pattern [20| and to the motion of remnant itinerant 
electrons not gapped by the ordering transition [16, 18]. 

A richer behavior is revealed by the data at finite tem- 
peratures, that clearly indicate the emergence of a tem- 
perature scale T* that marks an analogous suppression 
of the kinetic energy occurring within the homogeneous 
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FIG. 2: (color online) . Kinetic energy of the A^s = 12 interact- 
ing system at U /t = 5 (a) and U /t = 15 (b); T-dependence of 
low frequency optical weight (c), and integrated optical spec- 
tral weight I{uj) for T/t = 0.01, 0.03, 0.05 and 0.1 - 0.5 in 0.1 
intervals (d). Arrows in (a),(b),(c) correspond to inflection 
points of the curves, defining the temperature scale T* . 



metallic phase (we define T* as the locus of steepest vari- 
ation of K with temperature, denoted by arrows in Fig. 
[2|) . The scale T* appears to be entirely controlled by the 
approach to the zero-temperature ordering transition, a 
behavior that is strongly reminiscent of what is expected 
close to a QCP. The situation is similar at U /t = 5 and 
U/t = 15 [Fig. IH^b)], although in the latter case the 
kinetic energy ratio at T = is already reduced down 
to values K/Kq ^ 0.7 before entering the CO phase at 
Vco = 4.83i, which is indicative of a moderately corre- 
lated electron liquid. In this case the quantum critical 
behavior adds up to the correlated electron picture, af- 
fecting the motion of electrons that have already been 
slowed down by local electronic correlations. 

The above observations can be directly related to the 
low-energy quasiparticle properties by analyzing the tem- 
perature dependence of the integrated optical weight 
I{uj) = L a{uj')dijj' . The low-frequency integral I{ijJ = 
0.5i), reported in Fig. HJc), comprises most of the quasi- 
particles contributing to the Drude behavior in a nor- 
mal Fermi liquid, while excluding higher-energy incoher- 
ent excitations arising from the strong electronic inter- 
actions. Comparison of Figs. WL^) and (c) demonstrates 
that the strong reduction of kinetic energy above T* pri- 
marily originates from a drastic suppression of the low- 
energy coherent quasiparticles. Finally, Fig. [2][d) illus- 
trates I{ijj) at a given V/t = 4.8 just below the CO tran- 
sition, showing that the quasiparticle weight lost at T* 
is partly transferred to high-energy excitations, that are 
broadly distributed on the scale oiU.V. 

The above results are summarized in the finite tem- 
perature phase diagram of Fig. [3J that constitutes the 
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FIG. 3: (color online). Finite temperature phase diagram 
illustrating the emergence of non-Fermi liquid properties in- 
duced by the proximity to the charge ordering transition. The 
interaction strength has been rescaled to the critical values to 
compare QCP related properties for weakly and moderately 
correlated systems. 



central result of this work. Scaled units V/Vco are used 
so that the weakly {U/t = 5) and moderately {U/t = 15) 
correlated cases can be directly compared, illustrating 
the universality of the T* phenomenon in proximity to 
the CO instability. Both T* and Tco appear to vanish at 
Vco leading to a funnel-type 'bad' metallic region with 
strong quantum critical fiuctuations. 

To obtain further insight into the behavior of quasi- 
particles near the CO instability and make contact with 
the established concepts of quantum criticality, we have 
calculated the specific heat coefficient 7 — Cy/T. Our 
results, reported in Fig. HJa), resemble the behavior 
of nearly two-dimensional antiferromagnetic metals in 
which a singular increase is expected upon lowering the 
temperature close to the QCP, crossing over to a constant 
value at the onset of Fermi liquid behavior [J, 12^. Curves 
similar to those in Fig. ^a) are commonly observed in 



heavy fcrmion systems [J, |23|, \24\. As a striking confir- 
mation of the QCP scenario emerging from the preceding 
paragraphs, we see that there is a direct correspondence 
between thermodynamic and dynamical properties: the 
peak position in the specific heat essentially coincides 
with the temperature T* derived from the kinetic energy 
and the low-frequency optical integral near to the QCP 
(see Fig. [3]). 

Finally, we discuss the behavior of the effective mass 
as extracted from the low-temperature limit of the spe- 
cific heat coefficient (in practice we estimate m* from the 
peak value of 7 to overcome the numerical limitations of 
the Lanczos technique at low T). It can be expected 
on general grounds that the strong electron-electron in- 
teractions responsible for the 3-fold CO will strongly af- 
fect those parts of the Fermi surface that are connected 
by momenta closest to the ordering wavevector Q. This 




FIG. 4: (color online) (a) Specific heat coefficient 7 — Cv/T 
in units of 1/t for U/t = 15. The curves are shown down 
to a temperature comprising ^ 10 excitations in the Lanczos 
diagonalization. (b) effective mass ratio, m* /mi,, estimated 
from the peak value of 7 shown in (a) . 



should lead to the emergence of "hot spots" with diver- 
gent effective mass, m*/mb ex ln(l/|y — Vco\), in full 
analogy to the situation encountered in metals close to a 
magnetic instability [22, l25| . The effective mass reported 
in Fig. mjb) indeed shows a marked enhancement at the 
approach of Vco: that adds to a moderate renomalization 
m* /rai, < 2 provided by non-critical electronic correla- 
tions. Whether the destruction of quasiparticles remains 
confined to such 'hot spots' or spreads over the whole 
Fermi surface is an unsettled issue that is also actively 
debated in the context of heavy fermion materials [2, |3J . 
Concluding remarks. Our results indicate the occur- 
rence of non- Fermi liquid behavior driven by a combina- 
tion of electronic correlations and quantum critical fluc- 
tuations close to a CO instability in quarter-filled organic 
conductors. Quantum critical behavior has been recently 
reported in transport studies of the quarter-filled com- 
pounds K-(DHDA-TTP)2SbF6 and (MeDH-TTP)2AsF6, 
by tuning the system across the CO transition via an 
applied pressure [a, l7[. A stringent verification of our 
theoretical predictions could be achieved in the conduc- 
tor 61— (BEDT-TTF)l3, whose band structure is properly 
described by the model Eq. (1). Remarkably, this is 
the only salt of the 9 family exhibiting superconductivity 
and is at the edge of CO. Recent optical studies [26| have 
shown a rapid loss of electron coherence upon increasing 
the temperature, associated with a marked reduction of 
kinetic energy as obtained here, and the evolution of the 
integrated spectral weight /(w) with temperature found 
experimentally compares very well with our result in Fig. 
WiA) ■ The presence of an unexplained far infrared absorp- 
tion peak whose position is controlled by the temperature 
scale T alone could be a clue of an emergent collective 
excitation of the QCP 27|. This material undergoes a 
CO transition under pressure [28|, which could be di- 
rectly exploited to probe the quantum critical behavior, 
applying the plethora of experimental techniques that are 
commonly used in the study of heavy fermion materials. 
Hall coefficient as well as de Haas-van Alphen experi- 



ments appear to be ideal probes to test whether quasi- 
particles are destroyed over the whole Fermi surface or 
on some regions only, shedding light on the nature of the 
charge order QCP. 
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